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j We study the ionic equilibria and interactions of neutral semi-permeable spherical shells immersed in electrolyte solutions, 
including polyions. Although the shells are uncharged, only one type of ions of the electrolyte can permeate them, thus leading to 
{S| "a steric charge separation in the system. This gives rise to a charge accumulation inside the shell and a build up of concentration- 

, i dependent shell potential, which converts into a disjoining pressure between the neighboring shells. These are quantified by using 

^ -the Poisson-Boltzmann and integral equations theory. In particular, we show that in case of low valency electrolytes, interactions 
between shells are repulsive and can be sufficiently strong to stabilize the shell dispersion. In contrast, the charge correlation 
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effects in solutions of polyvalent ions result in attractions between the shells, with can lead to their aggregation. 
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In recent years, there has been much interest in creating self- 
assembled micro- and nanocontainers for various technolog- 
ical applications, including pharmaceutics, cosmetics, food, 
chemical and biotechnologies. In contrast to typical macro- 
scopic containers, the shells of microcapsules are molecularly 
thin and, therefore, it is usually hard if not impossible to make 
them completely impermeable for the solvent, the small ions, 
or low molecular weight solutes. Typical examples of such 
1 shells are vesicles and liposomes with ionic channels, 1,2 var- 
ious types of micro- and nanocapsules,-"— cell-^ and bacte- 
rial 9 " 11 membranes, viral capsids. 12,13 Such shells inacontact 
with a polyelectrolyte solution have been employed for mea- 
suring elasti c V 4 i 15 and elasto-plastic 16 modulus of the shell's 
material and for selective encapsulation.—"— Recent stud- 
ies have shown that mechanical properties of semi-permeable 
capsules filled with polyelectrolytes, their aggregation and and 
sedimentation stability are strongly influenced by ionic Don- 
nan equilibria between the container interior and the surround- 
ing solution.^ 

During the last few years several theoretical and simula- 
tion papers have been concerned with the Donnan equilibria 
in charged systems with semi-permeable walls. A large frac- 
tion of these deal with the calculation of ion profiles and the 
excess osmotic pressure on the membrane, by using the lin- 
earized mean-field Poisson-Boltzmann (PB) theory^"— and 
molecular dynamics simulations. 22 - 24 Other authors used in- 
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tegral equation theories to address the effect of the wall thick- 
ness on the electrolyte distribution. 25 Moreover, we have re- 
cently discussed a related situation, when two ionic electrolyte 
solutions consisting of large and small ions are in equilibrium 
with a thin film bounded by the flat semi-permeable mem- 
branes.— As counterions could escape from electrolyte solu- 
tions, at some separation their clouds overlap and, as we have 
shown by the non-linear PB theory and molecular dynamic 
simulations, give rise to a repulsive force between the mem- 
branes. However, many aspects of interactions of membranes 
have been given insufficient attention. 

In the present paper, we are considering a particular case of 
spherical semi-permeable shells in osmotic equilibrium with 
outer electrolyte solutions, including solutions of multivalent 
ions. We will use the non-linear PB equation and integral 
equation theory based on the Ornstein-Zernike (OZ) equation 
with the hypernetted chain (HNC) closure to evaluate the ionic 
distributions, accumulated charge inside the shell, and interac- 
tion between the shells. 



2 Model and Theory 



We consider semi-permeable spherical shells of radius R im- 
mersed in an electrolyte solution, as shown in Fig. Q] The 
electrolyte cations are characterized by charge Z + e, where e 
is the electron charge, and bulk concentration C+, while the 
anions by Z e and C~, respectively. The shell is uncharged 
and infinitesimally thin. For the sake of definition we assume 
that the shell is permeable for cations and impermeable for 
anions. Furthermore, we impose electroneutrality of the bulk 
electrolyte, so that Z + C+ = Z_C~. 
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Fig. 1 Schematic of a system consisting of two shells immersed in a 
polyelectrolyte solution. The shells (large spheres) are permeable 
for counterions (small spheres) only and impermeable for polyions 
(medium-sized spheres). 



2.1 Nonlinear PB theory for single capsule in an elec- 
trolyte 

We first describe an isolated neutral shell immersed in an elec- 
trolyte using the Poisson-Boltzmann equation. We introduce 
the dimensionless potential <p = y/e/ksT , where \j/ is the elec- 
trostatic potential, kg is the Boltzmann constant, and T the 
absolute temperature. The equations and the boundary condi- 
tions for the dimensionless potential are then given by 
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<p/|x=0 = 9i\x=l = 9o\x=\ 
<Po|*=«, = ^|^=i = ^|^ = i 



where x = 4, index i refers to the capsule interior, while o to 
the outer solution, K 2 = Anl B C+, ff„ = K^/Z + (Z_ +Z+), and 

the Bjerrum length is l B = AK £ ekgT ■ 

To solve the equation we map the far field point x = °° to a 
finite x* and then, with the substitution 



t,i = XiVx* - 1, § = Vx 
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we derive 
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= (x*-l)(fc/?) 2 Z + ( e z -'' 
with the boundary conditions set by 

<P/l|=O = 0, 9il5=V^ = ^l^=v^- 
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in the interval t, € [0; \Jx* — lj . 

At 4 = the first of the equations in (O is singular. To re- 
move the singularity we can use <p"|^= = - frarrft e~ z+<p 't Q ). 
Numerical solution of simultaneous equations (0 can be then 
found using the following procedure. First, as we map an in- 
finitely remote point onto a finite point x* , we find x* inf such 
that the potentials <p,- and <p„ as well as their derivatives remain 
constant on the interval % e [0; yx* — lj . The correspond- 
ing equations ^ with boundary conditions Eq. (01 can then 
be solved using the three-stage Lobatto Ilia formula (bvp4c 
function in Matlab). 

In the limit of large radius or high electrolyte concentra- 
tions, kR ^> 1, the shell is equivalent to a planar isolated semi- 
permeable membrane. The shell potential (p s is then defined by 
the boundary conditions (p'\ x =±oo = and takes the form 22 ' 26 
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Thus, in monovalent electrolytes, the induced Donnan poten- 
tial on the shell does not exceed <p* = ln2 s» 0.693. In a diva- 
lent electrolyte, its limiting value is <p* = jln3 « 0.549, and 
in trivalent electrolyte <p* = | In 4 m 0.462. At room temper- 
ature these correspond to 18 mV, 14 mV, and 12 mV, respec- 
tively. Note that while the PB theory predicts the vanishing 
potential for larger polyion valencies, the mean-field model it- 
self becomes increasingly inaccurate for multivalent ions. As 
the correlations come into play, the amount of charge the shell 
can accommodate will be even higher. In charged colloidal 
systems this leads to charge inversion, which, as we will see, 
appears also in our system. The charge inversion phenomenon 
is related to the charge discreteness and correlations, so that it 
lies beyond the mean-field description. 27 > 28 

We can calculate the accumulated charge inside the shell 

as Qn = Q/e = (jjjj <Po \x=i- A (neutral) shell with the accu- 
mulated inner charge will be referred below to as a (charged) 
capsule. We can now analyze an asymptotic behavior of the 
external potential <p„|. r =i at kR — > °°. It is easy to see that 



d<p 
dx 



{KR) 2 F[<p s }+4{KR) y/F\v\d<p (6) 



where we introduced a function F [<p] = 
2(f^[exp(Z_<p)-l] + [exp(-Z+<p)-l]). At KR » 1 
we have <p s w <p*, then in <p* we can omit terms of the order 
of (kR)~ . Similarly, we can omit the second term in Eq. © 
to find 

&H = rt-(r^VF\ri\ ( 7 ) 
In a symmetric electrolyte, Z — Z + = Z , Eq. (0 takes a sim- 
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Fig. 2 Geometry of a system consisting of two interacting capsules. 
The polar angle 6 and azimuthal angle <j) in the spherical polar 
coordinates, the center-to-center distance d and surface-to-surface 
separation h are shown. 
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which for a monovalent electrolyte reduces to = kR = 
KgR/y/2. Thus, for large capsules or at high electrolyte con- 
centrations the ion charge accumulated inside the shell is 
growing proportionally to the square root of the salt concen- 
tration (ff) and/or to the square of the radius. 

2.2 Interaction of capsules in the superposition approxi- 
mation 

In this section, we evaluate the interaction of two semi- 
permeable shells (capsules) in the superposition approxima- 
tion where an analytical solution can be found. 

We consider two spherical capsules of radius R with their 
centers located at (0,0,0), (0,0, d), as shown in Fig. [2] where 
d is then the distance between the centers such that d ^ 2R. In 
the first-order approximation, the resulting potential will be a 
sum of the individual potentials of the shells 22 



(p {x,y,z) = <pi X) (x,y,z) + q>o> (x,y,z - d) 
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This is justified provided the Donnan potential on the shell is 
small (<p — < lj. This assumption is reasonable for a 
long-range part of the solution. In other words, our arguments 
are valid for a relatively large distance between the shells (see 
AppendixlAlfor more details). 

For the shell centered at the origin, its surface potential can 



be rewritten in spherical polar coordinates as 



tf(e,d) = <p s 



0(0, d) 



exp(jQ,/f(l-#(M))) (10) 



where #(0 ,d) = \J 1 + (^) — 2 (^) cos and the surface po- 
tential <p s is calculated numerically as described in Section 2.1. 
Then, the pressure exerted by the non-permeating ions on the 
shell will take the form 



p(8,d) = ^k B TC+exp(Z-<p s (9,d)) 



(11) 



The z-component of the net force acting on the shell at the 
origin is given by 

{F e ff)z = Jj p(n l -n 2 )da = - JJ p-^da (12) 



where ni = — (x,y,z) /R is the unit normal vector to the cap- 
sule shell directed toward its center at (0,0,0), 112 = (0,0, 1) 
the unit vector pointing from the center of capsule at (0,0,0) 
toward the center of the second capsule at (0,0, d). From Eqs. 
( [Tol l and ( fTTT i we finally get 
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= InR 1 J p{B,d) sin0 cos Odd (13) 



It is now instructive to evaluate the asymptotic behavior of 
the repulsive force F e ff (r) for monovalent electrolyte in the 
long-distance limit, r^> R. This can be done by expanding the 
exponential function in the integrand 



o Z -<Po = 1 + ffl e 

r/R 
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where we replaced d with r for consistency with the rest of the 
paper. This yields for monovalent electrolyte at K„R 3> 1 



,(r) = 



InCty^R^e-^- 2 ^ ln(2)K fl 2 e -M'-2fi) 



K„r 21b 
An integration of the force Eq.(fT4b gives 

U etf (r) 27Cq> s e^C+R 2 e 2K ° H 



(15) 



k B T Ko 

\n(2)K t ,R 2 e 2K ° R 
2F B 



r(0,Kor) 
r(0,JC o r), (16) 



where T(0, K u r) is the incomplete gamma function. In the fol- 
lowing, we calculated U e ± ± (r) numerically using this equation 
and integration of the mean-field force [Eq. dl3H. 
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We remark that our asymptotic results differ from what is 
expected for (impermeable) charged colloids. Indeed, we have 
obtained a Yukawa-like long-range behavior for the force (not 
for the potential), which indicates the regime of an "ideal gas" 
pressure of overlapping ionic atmospheres. 29 Such a decay of 
the interaction force [Eq. (JT3J)] is likely a consequence of the 
two simplification used (a neglected effect of the second cap- 
sule on the charging and an assumption that kR^> 1). Nev- 
ertheless, as we will see below, this approach leads to quite 
accurate results at large separations. 

2.3 Integral equation approach 

In many experimental systems of interest multivalent ions are 
present. In this situation, effects of charge correlations on in- 
teractions between shells might be important. To address this 
issue, we need a more accurate treatment of the ionic distri- 
butions. Here we use the integral equation theory based on 
the OZ equation. To describe a capsule in an electrolyte, we 
will treat it as a separate type of particle, whose concentration 
is infinitesimally small. More generally, an inhomogeneous 
external field can be introduced as a new particle type at an 
infinite dilution. 25 i 3Q i 31 

A multicomponent OZ equation for n + 1 particle types has 
the following form 

n+l p 

M r 2i) = c, ; (r 2 i) + 9m I h im (r 2 3)c m j(r 13 )dr 3 (17) 

m=l J 

where p m is the number density of particles of type to, 
hij(v2i) — g//(r 2 i) ~ 1 and cy(r 2 i) the full and direct corre- 
lation functions for particle types i and j taken at r 2 and rj, 
gi;'( r 2i) is their pair distribution function, and r 2 i = r 2 — rj. 
The problem of calculation of the distribution functions from 
the known pair potential requires an additional relation for 
each pair ij, relating fty(r 2 i) and cy(r 2 i). Here we will use 
the hypernetted-chain closure, which is commonly used for 
electrolyte systems.— 

Cij(r 2 \) = -Puijfai) + fcyfoi) - ln £u( r 2i ) (18) 

where K,;(r 2 i ) is the interaction potential for particles of types 
i and j and J3 = 1 /kgT . 

For an n-component fluid in an inhomogeneous field, the 
field can be considered as the (n + l)-th particle type (we 
denote it by index 7), whose density is infinitesimally small 
Py — > 0. The full correlation function for the field, type 7, and 
particle of type j is given by 

n a 

hyj (r 2 1 ) = cyj (r 2 1 ) + £ p m I h y, n (r 23 ) c mj (r 1 3 ) dr 3 (19) 

m=l J 

The full correlation functions of the remaining particle types 
satisfy n-component OZ equations analogous to Eq. (fTTT i but 



without 7 and alow us to calculate the direct correlation func- 
tions c m j(rri). In this scheme, the pair correlation function 
gyj( r 2l) is equivalent to the inhomogeneous one-particle dis- 
tribution function gj(r\) of particle type j in an external field 
imposed by the capsule . Then, functions hyj (r 2 i ) and eyy (r 2 i ) 
can be replaced by hj(r\) — gj(ri) — 1 and c/(ri), respec- 
tively. Therefore, the local concentration of particles of type j 
in the external field becomes 



P;(n) =Pj8j( r h 



(20) 



Now we can use the HNC closure ( TT81 to substitute for 
c w -(r 2 i) inEq. (O 



g;(ri) = exp|-)3«y(n) + £ p,„ J h m (r 3 )c mj (ru)dr 3 j 

(21) 

where index 7 is omitted for consistency with Eq. (1201) . Here, 
the function 0,(^13) in the integral, Eq.tBTTi is replaced by 
the direct correlation function for n-component homogenous 
fluid, which is given by Eq. (TTTb without particles of type 7 
via the HNC closure ([T§i 

By solving Eq. (I2TI 1 together with the closure relation (1181 , 
we can find gj(r{) and c,-(ri). Hence, we can also calculate 
the distribution of particles j near the shell, pjif) — pjgj(r). 

Bulk concentration of particles of type j is set by ciP — 
lim Pjgj(r). Since gj(r) — >• 1 at r — >• °o, we have c£ ; ' = pj. 
Similarly to Eq.ffJTJ we can evaluate the pair distribution func- 
tion for capsules gyy(r 2 i ) in the limit p r — >• from 

n ,. 

]ngyy(rn) = -j3M yy (ri 2 )+ p m / n m (r 2 3)c m (ri 3 )Jr 3 

m=l J 

(22) 

where Uyy is the direct interaction potential for two capsules 
with their centers at r\ and r 2 . To calculate the potential of 
mean force for the capsules, we can use a one-component OZ 
equation (OCM) for a system containing capsules only such 
that it produces the same full pair correlation function as the 
full system, h OCM (ri 2 ) = hyy (ri 2 ). In the limit Py — > 0, we 
find that c 0CM (r 12 ) = h OCM (rj 2 ). The HNC closure d then 
gives pU e ff(r) = — Ingyy(r), which is the required potential 
of mean force. 

In the full system, we consider a two-component primitive 
electrolyte: positively charged spheres, which penetrate freely 
through the shell and mono or trivalent anions, which cannot 
penetrate into the capsule interior. The ion charge Z±e is put 
at the ion center. The interaction between cations as well as 
between cations and anions is set via a hard sphere potential. 
The solvent is considered as a continuous uniform dielectric 
medium characterized by the Bjerrum length lg =2. Ions of 
types i and j separated by a distance r between their centers 
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interact via 



u ij( r ) = 



l B {k B T 



QiQj 



r > a; 



(23) 



and dij for monovalent ions are a ++ = 



where i,j = +, — 

a + _ = a = a, while for trivalent ions a ++ = a = a, 

a + - = 2cr. Note that at low electrolyte concentrations the 
ionic distribution is sensitive only to a+_, which determines 
the cation-anion closest approach distance and, hence, the 
binding energy. The two other diameters are unimportant 
since like-charged ions rarely collide at low concentrations 
due to Coulomb repulsion. As for the cation-anion distance, 
we have increased the contact distance to a + _ = lg to reduce 
the ion association, and make the correlation effects more pro- 
nounced. Global electroneutrality was provided by setting 
Z + C+ = Z C . The interaction of anions with the capsule 
shell was described by the hard sphere potential 



— nonlinear PB 
— linear PB 




o. 



r<R 
r> R 



(24) 



We will further express all the length variables in a. 

We solved the multicomponent OZ equation with the HNC 
closure using the Ng's method with the modified version of 
program PLOZ. 32 



Fig. 3 Dimensionless electrostatic potential in the center of a 
capsule, (p m , for electrolytes with Z + = 1 and Z = 1 obtained 
within the nonlinear PB equations (solid curve) and the linearized 
PB equations 22 (dashed curve). We note that the curves for divalent 
and trivalent electrolyte, Z =2,3 would coincide with the result 
for monovalent salt on this scale and are omitted here. 



3 Results and Discussion 

3.1 Distribution of charges and potentials near a single 
shell 

Now we focus on a distribution of ions and of electrostatic 
potentials around a single shell. First we analyze the predic- 
tions of the non-linear PB theory and its linearized version. 
Then we compare the results of the PB theory with calcula- 
tions done using the integral equation approach. 

The distribution of the electrostatic potentials at the cen- 
ter of the capsule, <p,„, versus kR is shown in Fig. [3] One 
can see that the potential at the center grows rapidly with kR 
and reaches very large values above IOIcbT (or 250 mV) at 
kR ^> 1 (no overlap of the inner ionic layers). The computed 
data show no dependence of <p„, on the valency of large ions, 
Z . Note that the linear theory predicts a saturation of the 
potential at <p m = 1 (see 22 ) and, hence, significantly underes- 
timates the potential at the center at large kR. At kR <C 1 
(strong overlap of the inner double layers), the potential at the 
center vanishes. Such a situation would be realistic for very 
dilute solutions and/or very small radius of the shell. These 
results are qualitatively similar to those obtained previously 
for a thin gap between two semipermeable membranes.— 



Fig. [4] includes theoretical curves for the potential of the 
shell, <p s , calculated for several Z_/Z + , by using the PB ap- 
proach. At small kR all curves coincide, but at large kR the 
surface potential becomes sensitive to Z_/Z+ and decreases 
with this ratio, which is in agreement with the asymptotic re- 
sult [Eq.(f5]>] . Also included in Fig.|4]are the theoretical curves 
predicted within the linearized PB theory. At large kR there 
is a discrepancy between predictions of nonlinear and linear 
theory, which becomes smaller at larger Z_/Z+. An explana- 
tion for for these deviations is the different (underestimated) 
asymptotic value of the bulk Donnan potential in the linear 
approach. 26 

The build up of surface and center potentials is caused by 
a steric charge separation in our system. This is further illus- 
trated by the data in Fig. [5] where the net accumulated charge 
inside the capsule is shown. The charge grows linearly with 
kR when the latter is large, but is proportional to [kR) 2 when 
it is small. Note that the accumulated charge is not too sensi- 
tive to the valency of ions. The asymptotic law for the charge 
[Eq. (0] is also included (dotted-dashed curve), and becomes 
fairly accurate at large kR. 

Figure|6]illustrates the ionic distributions, g_ (r) and g+(r), 
in the vicinity of the shell, calculated for several concentra- 
tions of trivalent electrolyte by using the integral equations 
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Fig. 4 Dimensionless shell potential, (p s (KR), for electrolytes with 
Z + = 1 and Z_ = 1,2,3 obtained from the nonlinear PB equations 
(thick curves) and the linearized PB equations (thin curves). The 
limiting Donnan potential of an isolated shell [Eq.l(5] is shown by 
the dotted lines. 
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Fig. 5 Net charge inside the semi-permeable shell from the 
non-linear PB equation (thick curves), linearized PB equation (thin 
curves) and the asymptotic behavior (dotted curves) as predicted by 
Eq. 10 for an electrolyte with Z+ = 1 and Z_ = 1,3. 
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Fig. 6 Concentration profiles for (a) cations (g+{r)), (b) anions 
Q>- ( r ))> obtained from the OZ-HNC calculations (thick curves) and 
PB equation (thin curves) for an electrolyte with Z+ = 1 and Z — 3. 
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technique. We see that at small concentrations the PB equa- 
tion is very accurate, with only very small quantitative devi- 
ations from OZ-HNC. However, at large concentrations the 
agreement between the PB and OZ-HNZ is good only at very 
large r, but there are some qualitative discrepancies near the 
shell, both for anions and cations. The nonlinear PB approach 
leads to larger values of concentrations of multivalent anions 
at the shell than those predicted by OZ-HNC. Integral equa- 
tions theory also shows oscillations of the local concentration 
of trivalent anions and even counterions near the shell, which 
are an indication to ionic correlations in the system neglected 
in the mean-field theory. One can suggest that the layer of 
polyanions at the shell formed due to their attraction to in- 
ner counterions is charged stronger than the interior, so that it 
condenses some outer cations. The PB theory fails to detect 
this phenomenon, which is similar to the effect of charge in- 
version for charged colloids. 27 i 28 The criterion for the charge 
inversion and like-charge attraction is normally given by the 
value of the coupling parameter E = 2%a s l\Z^ > 10, where 
(7 S is the surface charge density.— It is satisfied at high va- 
lencies and high surface charge densities. This approach can 
be used to quantify roughly the phenomenon we observe, but 
we suggest to use the inner charge, (7 V = Q„/(4nR 2 ), to eval- 
uate the onset of ion correlation effects. Then the S values for 
shells of radius R = 200 and trivalent ions at the given range 
of concentrations varies from S as 1 at the smallest to E a 23 
at the largest concentration, so that it is not surprising that we 
have observed the correlation effects in the latter case. 



3.2 Interaction between capsules 

Fig. |7] shows the force acting to the shell as a function of inter- 
shell separation, h = r — 2R, calculated with OZ-HNC the- 
ory. The calculations are made for a monovalent electrolyte 
with fixed concentration by using several values of R. Also 
included are the predictions the mean-field [Eq. ( fT~3T >1 and 
asymptotic [Eq. ( TT3T )1 theory. The capsule interaction is re- 
pulsive and its amplitude grows with K a R, in accordance with 
predictions of Eq. ( [Tol l. At large K a R the mean field theory vis- 
ibly underestimates the value of a repulsive force. The striking 
result of our calculations is that approximate asymptotic solu- 
tion is surprisingly accurate, even at very short distances. 

The corresponding potential of mean force for capsules, 
U e ff(r), is presented in Fig. [8] Note that it reaches high val- 
ues at a contact, especially for larger shell radii. For typical 
colloid systems these values would be sufficient to safely sta- 
bilize the dispersion of solid particles. Since the radii of shells 
correspond roughly to 20-150 nm, we conclude that suspen- 
sions of the shells of radii larger than ca. 100 nm should be 
stable against aggregation due to an accumulated inner charge. 
Fig. [8] also includes calculations made by using an asymptotic 
formula, Eq.dToTl. It can be again seen that it is surprisingly 
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Fig. 7 Forces between capsules of radii R = 50, 200, 400, immersed 
in a monovalent electrolyte of concentration C+ = 10~ 3 . The 
corresponding screening parameters are KqR ~ 11,45,90. Thick 
curves are OZ-HNC calculations, thin curves show the predictions 
of the mean-field theory [Eq. J13H, and the asymptotic results [Eq. 
d 1 5H are shown by the dotted curves. 



accurate and can be used to evaluate the interaction potential 
in many important situations. 

The picture changes drastically in a solution of trivalent ions 
as it is seen in Fig. [9] At small electrolyte concentration the 
interaction is similar to that in monovalent electrolyte. How- 
ever, at large concentrations the interaction changes qualita- 
tively: an attractive force appears at intermediate distances. 
The criterion for like-charge attraction is the same as for the 
apparent shell charge inversion. Note that the range of the at- 
traction is the same for all capsule sizes and reflects the mean 
distance between the polyions within the "condensed" layer. 28 
We also remark that the attraction becomes stronger for larger 
capsules. 

The phenomenon of like charge attraction of semiperme- 
able shells is similar to that observed earlier for charged col- 
loids. It had been reported originally 33 ' 34 that the attraction 
arises only at certain concentration of the multivalent salt. 
Later it has been shown that this effect is generic and observed 
in mixtures of charged colloids and polyelectrolytes. 27 i 35 i 36 
As we discussed above, the potential of the shell decreases 
with the polyion charge (see Eq. (|5]l). At the same time, the 
energy of electrostatic polyion binding at the surface, is grow- 
ing as Z_<p s « InZ . Thus, the polyions will be stronger at- 
tracted to the shell, so that the ionic correlation criteria are 
easier to satisfy at the higher polyion charge. This was indeed 
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Fig. 8 Pair interaction potential for capsules of radii R = 50, 200, 
400 in a monovalent electrolyte of concentration C+ = 10 -3 at 
Ib = 2 as obtained from the integral equation theory (thick curves) 
and the asymptotic formula dl6t (dotted curves). The corresponding 
screening parameters are K„R ~ 11,45,90. 



found in the recent simulation study.— Hence, we expect a 
very pronounced attraction between the shells immersed into 
polyelectrolyte solutions and the aggregation of the shells as a 
result. 




4 Conclusions 



We have presented the nonlinear PB theory, including numer- 
ical calculations and asymptotic analysis, and OZ-HNC data 
on interaction of semipermeable shells immersed in an elec- 
trolyte solution. We found that electrolyte-mediated interac- 
tion between the shells in low valency electrolytes is always 
repulsive and is sufficient to stabilize their dispersion. How- 
ever, it can become attractive due to ion correlation effects 
already for a trivalent electrolyte provided its concentration is 
large enough. 



Fig. 9 Pair interaction potential for two capsules of radii 
R = 50,200,400 in a trivalent electrolyte of concentration 

= 3 • 10~ 4 (thin curves) and 3 ■ 10~ 3 (thick curves) at lg = 2 as 
obtained from the integral equation theory. 
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A Appendix: Asymptotic form of the electro- 
static potential at large distances 

Here we show that an accurate solution can be found if we 
solve the nonlinear equation for a shell potential, but use its 
asymptotic form, which is similar to the solution of the lin- 
earized PB equation 22 



(pi 1] = (p s -exp(-K (r-R)). 
r 



(25) 



It is known that the far field part of the solution of the non- 
linear PB equation for spherical geometry can be presented in 
the same functional form. In Refs. 37 ' 38 , it has been shown that 
the potential for charged colloidal particles calculated from the 
nonlinear PB equation can be presented at large distances as 



<p(r) = - exp[-Ko(r-R)] 



(26) 



where the coefficient B is given by 



B — 4tanh 



1 



1 



2k„R 



-tanh 



16 K„R Jo 



F(f)- 



/ 2 (/)sinh(/) 



K,R - In 



tanh(//2) 
tanh(<p s /4) 



■df 



(27) 



h = sinh (y) tanh2 (y) - tanh2 ({) +2hi 



F(f) 



cosh(<p. s /2) cosh(//2) 
sinh 2 (<p s ./2) sinh 2 (J/2) ' ~ tanh (<p s /4) 



In 



tanh (f/2) 

tanh(<p s /4) 
(28) 

tanh(//4) (29) 



A similar approach for calculation of effective interactions 
between charged colloids has been also used in the limit 
kR — » oo. 39 This procedure is associated with colloidal charge 
renormalization, as it changes the prefactor of the screened 
Coulomb potential. In our system, however, the potential is 
small, so that the second and the third terms of Eq. d27l i can 
be safely neglected, so that B is very close to (p s . In other 
words, the use of the potential in the form Eq. dZSt is well 
justified. 
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